Two-timescale adiabatic expansion of a scalar field model 
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The analysis of gravitational wave data may require greater accuracy than is afforded by the adia- 
batic approximation to the trajectory of and field produced by a particle moving in curved spacetime. 
Higher accuracy is available with a two-timescale approach using as an expansion parameter the 
ratio of orbital time to radiation reaction time. To avoid apparent divergences at large distances, the 
details of the method are important, especially the choice of the foliation, the spacetime surfaces on 
which the orbital elements are taken to be constant. Here we apply the two-timescale approach to a 
simple linear model to demonstrate the details of the method. In particular we use it to show that a 
null foliation avoids large-distance divergences in the first-order post-adiabatic approximation, and 
we argue that this will be true more generally for a null foliation. 



I. INTRODUCTION 

A major target of gravitational wave detectors, such as LIGO, but especially LISA[1[, is radiation from extreme- 
mass-ratio inspirals (hereafter EMRIs), compact objects, approximated as point particles, orbiting in the field of much 
larger gravitating centers, typically Kerr black holes. In the zeroth approximation these particles move through the 
Kerr background on geodesies characterized by orbital constants E (energy), L (angular momentum), and K (Carter 
constant), and the gravitational radiation generated by the particle is found with perturbation theory to first order 
—i\ in the EMRI ratio, the ratio of the particle mass to the background mass. 

In reality, of course, the orbits cannot be geodesies. The loss of energy and angular momentum to radiation requires 
that the orbits evolve. The evolution can be understood either as a reaction to the radiation, or as the result of a self 
force, the force on the particle due to the spacetime perturbations created by the particle itself. In either viewpoint, 
the deviation from a geodesic is proportional to the EMRI mass ratio, and for most EMRI ratios means that the 
deviation is small enough to justify what is called the adiabatic approximation: the particle is considered to be moving, 
at any moment, on a geodesic with orbital "constants" E,L,K that are functions of time as the orbit evolves. More 
important, the same approximation is made for the gravitational waveform generated by the moving particle: the 
waveform is to be that from a geodesic orbit characterized by E, L, K, with the three constants varying in time. 
In this adiabatic viewpoint the understanding of the orbit and of the waves it generates is shifted to the question 
. of the evolution of the orbital constants. The earliest method for doing this was to balance the loss of E and L to 
the radiation energy and angular momentum going out to infinity and down the horizon 2J. This method, though 
physically appealing, is not well justified except for circular orbits. Furthermore, it cannot be used to find the evolution 
of the Carter constant. The alternative, a calculation of the self-force driving the evolution, suffers both from the 
need for regularization of the divergent perturbation fields at the particle, and the gauge aspects of the self force. 
Some consistency of the balance and the self force methods has been found, when time averaging of the self force is 
used[U, although in a restricted gauge. 

Many technical details remain unresolved. In particular, the self-force can appear to be entirely gauge dependent 
on a short timescale 4]. Indeed, along with the issue of gauge invariance, the most pressing problem in this area is 
the question of the timescale on which the adiabatic approximation can be made to work. For gravitational wave 
analysis, what is needed is a description that applies over a large change in the orbital parameters. An outline of such 
a method has been proposed based on considerations of second order particle perturbation theory 

In order to formulate a justifiable scheme for an adiabatic approximation, a two-timescale analysis will be needed 
in which the adiabatic approximation is only the first step in an expansion of fields in orders of the ratio of a slow 
to a fast timescale. The mathematical framework for doing this was already inherent in Ref. @, in which it was 
shown that a consistent two-timescale approach involves several difficulties. First, the post-adiabatic effect is at the 
same order as nonlinear effects, and the nonlinear gravitational perturbations can result in both infrared-type and 



'Electronic address: mino@tapir.caltech.edu 
^Electronic address: rprice@phys.utb.edu 



2 



ultraviolet-type divergences that are not easy to regularize. Second, there is a weak violation of the gauge condition 
for the linear perturbations due to radiation reaction. The complications of those issues can obscure technical issues 
inherent in the basic idea of the two-timescale method itself. In particular, if not carried out in the optimal way the 
two timescale method can lead to apparent divergences. This is motivation for what is done in the present paper, the 
application of the two-timescale approach to a well-understood model: a linear scalar field coupled to a point particle. 

We use Minkowski coordinates {t, x, y, z} and we denote the spacetime coordinates collectively by Greek indices 
as x a (a = t,x,y,z), and we use the Einstein summation rule with the flat metric, r\ a p = diag(— 1, 1, 1, 1). Where 
appropriate we switch to spherical spatial {r, 9, </>} coordinates. When we refer specifically to the coordinate location 
of the particle we use the tilde, as in x a . We denote the spatial coordinates by Roman indices as in x l (i = x, y, z). 
We adopt geometrized units in which c = G = 1, so that mass, length and time have the same dimension. 

As a simplest toy model for the practical demonstration in this paper, we calculate the scalar field $ induced by a 
point charge q according to the flat spacetime wave equation, 

□$(ar) = p{x) , p{x) = q [ dr<5 (4) (x - x{r)) , (1.1) 



where we use the box-operator as □ = — (dt) 2 + (d x ) 2 + (d y ) 2 + (<9 Z ) 2 . Here S^\x) is the 4-dimensional Dirac's delta 
and t is the proper time of the orbit. 

For simplicity, we assume that the particle motion follows a Newtonian second law for the potential of a central 
gravitating mass M, 

S W =~^- +a ' (L2) 

where v l — dx 1 /dt is the ordinary spatial velocity (coordinate change per coordinate time) and a 1 is the ordinary 
accerelation (in terms of coordinate time) due to scalar radiation reaction. 

One of the challenges of point particle perturbation theory is to extract a finite a 1 from the scalar field which is 
divergent along the orbit [?J • But our main interest in this paper is the calculation of the first post-adiabatic term in 
a two-timescale expansion. For that we shall not need an explicit form of a 1 , and will treat it as a known quantity. 
As one additional simplification, wc will consider in this paper only the case of quasi-circular orbits (orbits which are 
circular except for radiation reaction) . 

The structure of this paper is as follows. In Sec. UH we calculate the orbital equation including the radiation 
reaction effect. In Sec. IIII1 we give the adiabatic approximation for the scalar field, the leading order of the adiabatic 
expansion. Though the result here is familiar, the approach in this section clarifies the geometric meaning of that 
approximation and, in particular, introduces the idea of the spacetime foliation function for evolving orbital elements. 
In Sec. IIV1 we calculate the first post-adiabatic correction of the scalar field. We conclude and summarize in Sec. [Vj 
and discuss applications of this new expansion scheme to other problems in relativity. In order to have the main text 
focus on the central ideas, we have relegated many of the details to a set of appendices. 

II. ORBITAL EVOLUTION WITH RADIATION REACTION 

In the adiabatic expansion orbits are characterized by the orbital elements, the integral constants of the orbital 
equations in the absence of radiation reaction. The evolution of the orbit due to radiation reaction is then described 
by the evolution of the orbital elements^. This means that the orbital elements become functions of "time," i.e., of 
some parameter / along the orbit. It is through this dependence that the effects of radiation reaction will appear in 
the orbits. One of the central issues of this paper is to consider just what time slicing of spacetime is appropriate 
to the evolution of the orbital elements. That is, how is the orbital parameter / to be promoted to a function of 
spacetime location f(x a )7 

We align our coordinates so that the quasi-circular orbit is at z = (Cartesians) or 6 = it/ 2 (sphericals). In order 
to calculate the orbital evolution p.2| , including the radiation reaction effect, it is convenient to use the orbital 
energy e and the z-component of the angular momentum i defined by 

1 . . M 

e=-i)V , l = xv v -yv x . (2.1) 

2 f 

These are constants of motion which are conserved along the orbit in the absence of radiation reaction. Due to the 
radiation reaction, e and £ are not constant and the evolution equations are 

d dt dt j i d dt e dt „ 

e ~ a e = —v % a\ —£:= —a* = —(xa v - ya x ) . 2.2 
df df df df df df 
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Here we have computed the dynamics of the orbits with coordinate time t, but have used dt/df, evaluated at the 
orbit, to infer the influence of the dynamics on the orbital elements. We have also defined two components of 
radiation-reaction driven acceleration: a e and a 1 . 

Because radiation reaction effect is weak, we first integrate the orbital equation of motion in the absence of the 
radiation reaction. The orbital coordinates are written as 

r(f) = rc, 0(/)=tt/2, «£(/) = fi [t(f) - 1(/ = 0)] + <f>c , (2.3) 

where rc and f2 are the orbital principal elements, and <pc is the orbital positional elements. The function t{f) is 
obtained from the definition of the foliation function, / fi, re, 7r/2, <^>(/)J = f. Since the orbit (in the absence of 
radiation reaction) is circular, e and £ are related by 

M 



(2.4) 



It has been shown that this circularity relationship is not only correct in the absence of radiation reaction, but holds 
to all orders in an adiabatic expansion, as long as the radiation reaction is small 9]. From this it follows that the two 
components of the radiation reaction acceleration always obey 

<t = -li a °. (2.5) 

For the circular orbit all kinematical quantities are related. It is convenient to express them all in terms of the 
central mass M and the orbital speed v: 

1 2 n M M ^ V 3 . „, 

e = -^v\ i=~, r c = —, n = —. 2.6 

2 v v z M 

With these relationships we can regard v as an orbital principal element that governs all other principal elements. 
The evolution of v itself is 

d 1 di v 2 dt t 

— v = a = a . (2.7) 

df vdf Mdf V 1 

For the relationships in (|2.6|) to hold in the presence of radiation reaction, Q, must be defined to be Vl :— d<fr/dt. 
According to (|2.3p . this requires that 4>c evolves according to 

|<fc = -[i(/)-t(/ = 0)]|fi. (2.8) 

Because the evolution of the orbital elements is derived perturbatively after we calculate the scalar field by the 
adiabatic expansion, we start with the orbit as a function of the orbital elements, and we write the source term of the 
scalar field equation as p(x a ) — p (x a \C(f(x a ))) where we denote by C the collection of orbital elements. 

III. ADIABATIC EXPANSION - LEADING ORDER 

The adiabatic expansion can be elegantly formulated as a two-timescale expansion: One scale is the dynamical 
scale of the system denoted by Td yn , such as the orbital period or the wavelength of the scalar field. The other is the 
radiation reaction scale denoted by T ra( j, such as the timescale for a change in orbital energy due to radiation reaction. 
We assume that the radiation reaction scale is much longer than the dynamical scale and we use the dimensionless 
small value \i := Td yn /T ra( j < 1 as an expansion parameter. That is, we assume that the orbital principal elements, 
the elements directly related to v, evolve slowly on the timescale of the radiation reaction. To clarify the order of the 
two-timescale expansion, we can use p as the expansion index and can replace the /-derivatives of v by 

From (12.61). we have that 



*j ^"U> v (3 - 1} 



d d d 2 d 2 



rc^-Sl^O^), —rc~-^n~0(p 2 ), (3.2) 
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and, from (|278|) . that 



(3.3) 



One may ask whether the first /-derivative of the orbital positional element d4>c/df may become large for large t(f), 
invalidating the expansion. However, as we will see in the next section, this is not a problem. 

In the adiabatic expansion we regard the scalar field as a function of these orbital elements and we expand the field 

as 



$ (x\C(f(x<*))) = $<°> (x a \C(f(x a ))) + (x a \C(f(x a ))) + ^ (x a \C(f(x a ))) + ■■■■ 



(3.4) 



To see how this expansion is used, we consider the wave operator acting on the scalar field. When the differential 
operator acts on the orbital elements of the scalar field, it results in a factor fi according to the scheme of (|3.f I) . To 
the leading order (i.e., order /z°) of the adiabatic expansion, we have the equation for <I>( ) to be 



□$(°)(:c«|C) 



c=c(/o«)) 



= P (x a \C(f(x a ))) 



(3.5) 



On the left hand side of (|3 - 5(1 , the wave operator is considered to act only on the direct dependence on x a , not on the 
dependence implicit in the dependence of C on /. In this sense, then, the orbital elements can be considered to be 
constants. By also treating the C(f) on the right as if they were constants, we interpret the wave equation (|3.5[) as 



n$W(x\c) = P (x\c)) , 

for which we can give the well known solution for a circular orbit with constant orbital elements, 



(3.6) 



(x a \C(f(x a ))) = qyj 1 - r 2 c n 2 ^2(-imQ) ( h[ l) (mOr) j; {mQr c )0(r - re) 



+j l (mnr)h< t 1) (mnr c )6(r c - r))Y lm (e^)Y l * m (7r/2,0)e- m ^ t+ ^ , (3.7) 



where the orbital elements are evaluated by rc = rc(f(x a )), O = il(f(x a ))) and <fic = <l>c(f(% a )))- 

We note that, to this leading order of the adiabatic expansion, the foliation function does not play a crucial role in 
the practical calculation, but it does play an important role in the meaning of the expansion, in particular how the 
evolving elements of the orbit affect the solution at points off the orbit. It tells us that the solution at field point x a 
is the circular orbit solution for the values of rc, ft and <pc that occur at the "time" f{x a ) on the orbit. 



IV. FIRST POST-ADIABATIC FIELD 



In the field equation (|3.6[) . we ignored the derivative acting on the orbital elements and, by the rule (|3.ip . we found 
a field $ ( -°' ) that satisfies (jl.f [) only to accuracy 0(/i). We now turn to the calculation of the solution to the next order 
in /i. 

Using the orders for the derivatives of the orbital elements (|3.2p and (|3.3[) we can substitute the adiabatic expansion 
(|3.4p into the field equation and we can extract the terms of order fi. The details are given in Appendix A, and 

the result is a field equation for the first-order post-adiabatic field 



□$ (1) (x|C) 



C=C(f) 

P {1) (x\C(f)) 



(4.1) 



d a d a ^°\x\C) 

-d+a*® {x\C{f))g aP d a fd p f 



C=C(f) 



f P 3pf + da*™ ( X \C{f)) Df 





'dC a ~ 


} 


[ df \ 



(1) 



dp 



(1) 



(4.2) 



where the notation [dC/df]^ and [d 2 C / df 2 ] {1 \ means that these derivatives are only to be taken to first order in 
/i. On the left hand side of (|4.2p the derivatives in the □ operator are understood only to act on the explicit x a 
dependence, and not on the orbital elements. The first-order nature of the left hand side is due to the fact that 
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itself is defined to be first order in (|3.4[) . The result in (|4.2[) . with the orbital elements treated as effective constants 
on the left is the scalar version of the post-adiabatic equation first derived in Ref. [6j . 

The foliation function f(x a ) played only a passive role in the adiabatic approximation of Eqs. (|3.5|) - (|3.7p . By 
contrast, here the post-adiabatic terms depend explicitly on the choice of f(x a ). The choice of the foliation function 
is similar to a gauge choice, in that the value of any physical quantity cannot depend on the choice [lOj. As a useful 
check, the foliation-invariance of the first post-adiabatic approximation field $W is verified in Appendix iBl 

Though the choice of foliation function is arbitrary in principle, an inappropriate choice will make the first post- 
adiabatic term $W divergent so that a regularization procedure will be needed to extract physical information. In 
Appendix [C] we study the asymptotic behavior of the source for the first post-adiabatic term and we find that, if the 
foliation is asymptotically null, the first post-adiabatic term behaves as "J^ 1 -* — » 0(l/r) at large radius. 

Another insight on the appropriateness of a null foliation comes from the retarded time Lienard-Wiechert solution 
[H to (TO) : 



q_ 

47T 



(4.3) 



where £(r) has the components z a (r) of the orbit of the source particle, and where t' is the retarded time, a solution 
of 

t'=t-\x-at')\. (4.4) 
For the evolving equatorial orbital motion described by parameters ^(t), r c{ T ), and 4>c(t), eq. (|4.4[) becmes 

(4.5) 



t' = t- \Jr 2 + r%- 2rr c sin 9 cos (fl(t')t' + cf> c (t') 



The exact solution of (I4.5[) for t' as a function of t, r, 9, (j) is a null foliation, and would, in a sense, be the perfect 
foliation for the problem since with this foliation the solution (i.e. , that of (|4.3p ) is guaranteed to have no apparent 
singularities except at the source points. Equation (|4.5[) cannot in general be solved exactly, but for r » rp, and 
\t — r\ rc, the approximate solution is a simpler null foliation if = t — r. 
For simplicity, as well as from the form of the Lienard-Wiechert solution, we choose the foliation 



f(x a ) = t-r. 



(4.6) 



With this choice, □/ = —2/r in the second term of (|4.2| . g a/3 dafdpf = in the last term, the time coordinate of the 
orbit is given as t(f) = f + rc(f), and the evolution of cf>c is given by 



^c = -[t-r + r o (J)-ro(f = 0)] j f V. 



(4.7) 



This puts the source term of (14. 2[) in the form 



d t H — d r r 
r 



dC a 



$ (0) (x|C) 





'dC a ~ 


- C=C{f) 


[ df \ 


r — rc)-\ 





(1) 



Im K 



dv 

W 



(i) 



(4.8) 



where 



c(+) _ q(l)3 
J lm ~ J lm 



z—mVLr 



-s, 



(2)j 
Im 



s, 



(-) 



s, 



(l)h 



2hP(z) + 2z^-h^(z) - l 2zh ( , 1 \z) + l l l±H h t 1 \z) 
dz z 

3l{z) + -zTz Z3l{z) , 

z—m\ Ir 



(4.9) 



z—m^lr 
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c(p) 


2 








c(2)fc 


= —kAmv 

V 



2ji(z) + 2z—ji{z) - i2zji{z) + >-ji{z) 

az z 



Sim 1 — 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



with k representing either or j . 

For the most part, the derivation of (14. 8|) - fj4. 13[) is straightforward, but there is at least one point worth mentioning: 
The derivative with respect to 4>c entails a factor, shown in (|4.7[) . of t — r + rc{f) — fc(f = 0). The leading i will 
cancel with the £ arising from the derivative of e - imf2 (/)* j with respect to /, and the — r is rewritten as — mfiz, leaving 
us with 



rc(/)-r c (/=0) = n- 



1 - 



M/=o) 

ro(/) 



1 - 



^ 2 (/) 



« a (/ = 0) 



(4.14) 



Here the relationship tq — M/v 2 , from (|2.6p . has been used and, as usual, the symbols 17 and v, with no explicit 
argument are understood to be Cl(f) and v(f). 

For the source in (|4.8p - (|4.13p . the solution of (|4.1|) for <&W can be found using the standard retarded Green 
function for the wave operator. This solution is summarized in Appendix [D] It should be noted that the source 
term above has a discrete spectrum and contains only Fourier components with u> = mQ and that each (Zm)-modc 
of Sf m contributes only to the (im)-mode of <f>(°h This is due to our spherically symmetric choice of the foliation in 
(|4.6|) . With this choice, no explicit angular dependence is introduced into the differential operators. The separation 
of Fourier modes stops when we come to the explicit evaluation of [dv/df]^. This is a nonlinear step that mixes the 
contributions from all Fourier modes. 



V. SUMMARY AND DISCUSSION 



We have applied the two-timescale method to the simple well-understood problem of a slowly evolving quasi-circular 
orbit of a scalar charged particle in a central potential. The question motivating this work was that of the justification, 
and applicability of the adiabtic approximation and its extension. We have focused on the manner in which the slowly 
changing orbital elements of the trajectory, the radius ret angular velocity f2, and phase constant <pc are to be 
promoted to spacetime fields through the choice of some foliation function f(x a ). We have shown that this choice, 
which plays only a passive role in the lowest-order adiabatic expansion, is crucial for the post-adiabatic fields. In 
particular, we show that the first-order adiabatic field $W is well behaved as r — > oo if / is chosen to be t — r, but 
is not well behaved for a general choice of foliation function, in particular for the acausal choice / = t, which has the 
orbital elements change throughout spacetime according to coordinate time. 

An advantage of our simple model is that we have an alternative approach to the description of the field, the 
Lienard-Wiechert solution. We can get insight into the pathology of the wrong foliation by considering the r rc 
form of the Lienard-Wiechert solution in (|4.3[) 

$ = --<L_ 7 i . (5.1) 

A7Tr 1 - at') ■ x/r 

We have already argued, on the basis of the retardation equation (|4.4p or (|4.5p . that / = t — r is the correct large r, 
foliation function, but now we can consider what would happen for the foliation choice t' — t. In this case, at large 

r, the value of £(t' = t) would differ from the correct value £(t! = t — r), by an amount that is proportional to r and 
thus would constitute a contribution to the field that does not fall off as 1/r. 

It should be noted that the Lienard-Wiechert solution is itself a demonstration that for our simple problem the 
adiabatic approximation can be made to be exact over an arbitrarily long timescale. In (|4.3p one needs only to take 
the exact solution for t' as the foliation and to describe the motion £(r) in terms of evolving orbital constants, a 
description that can always be given if enough constants are included. From this point of view, adiabatic expansions 
of different orders represent various accuracies with which the foliation function / is equal to the t' that solves (I4.4[) . 
In our equatorial orbit problem, the usual lowest order adiabatic approximation (say that of (|3.7p ) corresponds to the 
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following procedure: One fixes the orbital constants tt, re and 4>C then solves symbolically for t'(x a ) parameterized 
by these constants and writes the Lienard-Wiechert solution in terms of that solution for t' . In this solution, one then 
promotes the "constants" to functions of whatever foliation is chosen. 

The Lienard-Wiechert solution, of course, applies only to problems involving linear fields obeying the flat spacetime 
wave equation. The two-timescale approach, in principle, applies much more generally. But at least four important 
questions remain about the two-timescale approach. The first is whether it converges for our simple problem (or 
similar linear problems). We feel that the existence and form of the Lienard-Wiechert solution all but guarantees that 
the post-adiabatic expansion, if done correctly, must converge. The caveat here, that the expansion must be done 
correctly, points to the need for the correct choice of foliation. A general argument based on (|5.1[) , and the details of 
the post-adiabatic equation in (|A9p and its higher order extensions, both point to the choice / = t — r being adequate 
to all orders. 

The second important question about the expansion is whether it converges for fields in curved spacetimes, such as 
black hole spacetimes. The issue is the choice of the foliation around the black hole horizon and the convergence of the 
post-adiabatic expansion. The curved-space generalization of the Lienard-Wiechert solution is obtained by the formal 
Hadamard series expansion using the bi-tensor formalism [T^. On the basis of this we suspect that the post-adiabatic 
expansion must converge with the choice of null foliation to the future horizon, but this is a question that deserves 
further study. 

The third important question about the expansion is the issue of the gauge condition when we apply this expansion 
to the gravitational perturbation problem[6(. In Sec. IIII1 we found that the first /-derivative of the positional 
element (13.31) could be large for large t(f). However, as we see from (|4.8[t . this does not invalidate the post-adiabatic 
expansion in our simple scalar model. This is not the case for gravitational perturbations and it is necessary to choose 
an appropriate gauge condition for the post-adiabatic expansion so that gravitational perturbations become valid on 
a longer time scale [6|. 

The last important question about the expansion is whether it converges for more realistic problems in relativity, 
especially nonlinear problems including the dynamics of geometry. Here we feel that the issue of divergences at large-r 
will probably be similar or identical to that for the linear problems, since the fields are weak at large r. The issue 
of the post-adiabatic expansion near strong field sources is one that will have to await the exploration of a nonlinear 
model. A point deserving particular attention is that the dynamical horizon may complicate the analysis and the 
convergence due to the time evolution of the black hole mass. To get some insights about nonlinear issues, we hope 
to apply the two-timescale method to is a post-helical solution based on a helical nonlinear problem, such as that in 
the post-Minkowski solution of Beetle et al. [l3[ 
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APPENDIX A: ADIABATIC EXPANSION OF THE WAVE OPERATOR 

In this appendix, we consider the two-timescale expansion of the wave operator. For the application of the two- 
timescale expansion, we regard the scalar field $ as a function of the spacetime coordinates {x a } and of the orbital 
elements, fi, rc and <j>c, collectively denoted C a . In addition to its direct dependence on the spacetime coordinates, 
the scalar field has an additional dependence through the foliation function f(x a ): 

$ = <&(x|C(/)) . (Al) 

As explained in (|3.ip , in the two-timescale expansion we assume that the derivative of the orbital principal elements 
is small, and in particular that 

Ml 

~ 0(f) , (A2) 

where fj, is the ratio of dynamic to radiation-reaction timescales. The derivative of the orbital positional elements is 
obtained from the orbita equation of motion. 

It is crucial to deal separately with the coordinate derivative of the scalar field with fixed C a and with the coordinate 
derivative of the scalar field due to the spacetime dependence implicit in C a (f(x a )). For clarity of the notation, we 
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[d^{x\C)] c=CU) , 

to denote the derivative with respect to the cordinates with the orbital elements C a fixed. We also use 

d 



dC a 



(A3) 



(A4) 



C=C(f) 



to denote the derivative with respect to the orbital elements C a . 

With these notations, the partial derivative of the scalar field is written as 

dC a 

d a $(x\c(f)) = [d a ^] c=c{f) + d a ^—d a f. 

The second partial derivative becomes 

dC a dC a 

d a di3$ (x\c(f)) = [d a d^] c=c{f) + [d a d a ^] c=cif) —d f + [d a df,*] 0=c(f) —d a f 



(A5) 



df 

dC a dC b 

+d a d b ^^-—d a fd f + a a $ 



df 

O a O(3f + —TprOafOpf 



df 



dp 



(A6) 



where we use the fact that the coordinate derivative with the fixed orbital elements commute with the derivative with 
respect to the orbital elements , i.e. , [dada^] c= c(f) = [®ct®a$\c=C(f)' 

Using the results above, we have that in general spacetime coordinates the wave operator acting on the scalar field 

is 



[□*] 



C=C(f) 



+g af3 (2 [d a d a ^] c=c[f) dfif + d a $d a d f - ri p d a ^f 

9 aP d a .fd p f, 



dC a 
~df 



aa dC a dC b a d 2 C a 



df df 



df 2 



(A7) 



where fTg is the Christoffel symbol. 



In applying (|A7[) to find adiabatic fields of any order, one must count the order by adding the index n 
of the field to the number of C a factors, keeping in mind that d 2 C a /df 2 can be first or second order. 

Thus, for example, g a/3 d a d a ^ 2 ^(dC a /df)f is third order and contributes only to the source for □<I>( 3 ), while 
d a d a & 2) {d 2 C a /df 2 )g a Pf, a f^ is a source for both and U^K 

To first order, the source that result from (|A7|) is that in (|4.2p . To second order, the result is 



n^ 2 \x\c) 



C=C(f) 

P {2) (x\C(f)) 



P {2) (x\C{f)) 



(A8) 



g«f>2 



dC a 



(i) 



df 

d a d a $M(x\C) 

d 2 C a 



dC b 



(i) 



8 a ^ 



d 2 c^ {2)S 



df 



g ap d a fd p f 



C=C(f) 



ddf + da*™ (x\C(f))Df 





~dC a ~ 


(1 


} 


[ df _ 





dp 



(1) 



g ap d a fd fi f. 



(A9) 



It should be noted that for null foliations the factor g a ® f a fp vanishes, simplifying the results. 



APPENDIX B: INVARIANCE OF THE FIRST POST-ADIABATIC EXPANSION 



Our adiabatic expansion uses a foliation function that is arbitrary in principle. If the result of the adiabatic 
expansion depends on the choice of this foliation, it means there must be a physical condition to determine the 



9 



foliation. Otherwise, the adiabatic expansion would not have physical predictability. In this appendix, we shall prove 
that the field found by the adiabatic expansion does not depend on the choice of the foliation to the first post-adiabatic 
expansion, and therefore, at least to this order, does have physical predictablity. 
We consider an infinitestimal change in the foliation function 



/(a*)- f(x a )+Sf(x a ), 



(Bl) 



and we show here that the first post-adiabatic field $W does not change to first order in Sf. 

There are two constraints on our choice of /, and hence on Sf. First, we recall the role of the foliation function. 
The orbital elements C, that are defined by the evolving orbit are promoted to spacetime functions by making them 
functions C(f(x a )). We choose to leave consant the value of / on every point of the orbit. (To do otherwise would 
not really involve a change in the foliation, but rather would be a change in the parameterization of the orbit that 
would change the functional forms of the C(/)). This gives us the constraint 



Sf(x a ) = , for x a on the orbit . 
Our second constraint is that we choose Sf to behave at infinity so that 



dSf_ 

dx a 



constant , as r — > oo . 



(B2) 



(B3) 



The reason for this, to be clarified below, is connected to the large-r boundary conditions that determine the solution 
for$W. 

The differential change Sf affects the field $ = + + • • • in two ways. First, at a given spacetime point x a 
the field formally denoted as $ = $' ' will have a change of order p according to 



$w (x\c\f)) -> *(°) (x\c(f + s.f)) = *(°) ( X \c(f)) + a$w (x\c(f)) , 



where 



Because dC a /df is 0(ji) by the ansatz of the two-timescale expansion (|A2 
post-adiabatic term. We note that 



&*<® (x\C(f)) = d a ^^Sf. 



(B4) 



(B5) 



□<5$ (0 >~ 


= □ 




C=C(f) 



dM^Sf 
df 



2g a ^d a d a ^S f + d a $wn6f 



§<&(°) is regarded as the part of the first 
dC a 



df 



(B6) 



The term d a U^ 5 f {dC a / df) has been omitted, since D^ ) vanishes except on the orbit, where, by constraint (|B2 
the factor Sf vanishes. 

To find the change in the formal post-adiabatic term $W we consider the transformation of p^ 



5p { 



i) 



(B7) 



We are interested in Sp^\ the change, to first order in p, to the right hand side of (|4.8|) induced by Sf. It is crucial to 
see that Sp^ does not come from changes in $^ ' nor to [dC a / df}^ in that right hand side; such changes would be 
post-adiabatic terms higher order in p. To first order in p the only changes induced by the change (|Blj) are those in 



SpM 



In principle, we can find t)^ 1 ) from 



d a d a ^°\x\c) 

'lC=C(f) 

+d Q $(°> (x\C(f))DSf 



g af3 d 5f 





'dC a ~ 


(0) 


} 


[ df _ 





c=c(f) 

Since the right hand sides of (|B8|) and (]F36|) add to zero, we have 



8p^ 



c=c(f) 



= 0. 



(B8) 



(B9) 



(BIO) 
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so the wave equation is sourceless for the change 5& of the scalar field, to first adiabatic order. 

In order to conclude that the change <5<I> is zero to first order, we need to rule out "free wave" solutions of (|B10I) . 
This can be accomplished by imposing the constraint in (|B3|1 . If some such constraint is not imposed then nonzero 
solutions for <5$ exist, and a change in the foliation function appears to make a change, to first adiabatic order, in 
the solution for But the change in the foliation function would entail a change from a solution that is causal, i.e., 
one that has only outgoing waves, to one that is not causal. In the adiabatic expansion we use only causal Green 
functions to solve for the various orders of <!>("). 



APPENDIX C: ASYMPTOTIC BEHAVIOR OF THE SOURCE TERM FOR THE FIRST 

POST-ADIABATIC TERM 



In this appendix, we consider a general spherical foliation defined as 

f = f(t,r), 



(CI) 



and we investigate the asymptotic behavior of the source term for the first post-adiabatic field with this general 
foliation function. We use the asymtotic form of the adiabatic scalar field written as 



* (0) (*\C(f)) = - J2 A lm ( C (f)) Y ?rn( 



)J_ e -im(n(/)(t-r)+<M/)) 



(C2) 



(See (I3.7P for the full expression.) 

The source term for the first post-adiabatic term, for the spherically symmetric foliation is 



c a =n,4> c 



2d t fd t - 20rfd r - \f + d 2 t .f - 8 2 r A ^$(°) (X. C ) 



dC a 



a-l(l) 



df 



+ ((d t f) 2 -(d t ff)^- a ^Hx,c) 



~d 2 C a ~ 




I df 2 _ 


1 



(C3) 



The regular behavior of $, and hence of the first post-adiabatic term <J>W is ~ 0(?* _1 ) at large r. For this to be 
the case, the first-order source term p^> must fall off at large r as O (f*~ 3 ) , and this condition constrains the foliation 
function. To order 1/r, in the expression above 



2d t fd t - 2d r fd r = -2im(d t f + d r f) . 



(C4) 



The function d a C a & ^ (x, C) falls off as 1/r, hence is asymptotically 1/r unless the leading terms vanish, that is 
unless 



= t f + d r f = d 2 f - 8 2 f , 



(C5) 



which indicates that the foliation surface must asymptotically be a future null cone. By these conditions, we have the 
foliation function written as 



f(t,r) -> f(t-r) at r->oo. 



If we put / = t — r in (|C3[) we get 



p«(x;C) = 2 J2 

c a =n^ c 



dC a 



4f 



(1) 



Since 



r r 
we have that p^'(x; C) falls off as r -3 , as it must. 



(C6) 



(C7) 



(C8) 
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APPENDIX D: INTEGRAL FORMULA FOR THE FIRST POST-ADIABATIC TERM 

We summarize the integral formula for the held equation under the outgoing-wave boundary condition. We consider 
the scalar wave equation 



aF(x a ) = S(x a ). 



(Dl) 



for flat spacetime, i.e., □ = —df + d 2 x + dy + <9 2 . In spherical coordinates, it is covenient to use the decomposition 



F{x) = Y,F lmu (r)Y lm (0, 



)(' . S(x) = Y,S lmu ,(r)Y lm (d,j>) t 



Imu Imuj 

A formal solution of the radial mode functions is obtained by the Green's method as 



Fi mu {r) = I dr'g lmu (r,r')S lmu (r), 

where 8(x) is the step function, i.e. 8(x) — 1 when x > and 9{x) — otherwise. 

For the calculation of the first post-adiabatic term, we write the source term schematically as 



(D2) 

(D3) 
(D4) 



Si mu (r) = s p lmu r c 5(r - r c ) 



Mi) 

Imuj 



1 zdz 



+ s 



h(2) 
Imuj 



, h(3) 



f4%) + z±h?\z) 



>9{r-r c ) 



Imuj 



ji(z) + -^-zji(z) 
z dz 



J (2) 
Imuj 



j l{z) + z -j l{z) 



n f \ 1(1+1) , , 

-i2zji(z)+i^ ! -ji{z) 



>6(r c - r) , 



(D5) 



where s 



h/j(l/2/3) 
Imuj 



does not depend on r. 



The radial mode function induced by the first term on the right hand side of (|D5|) is obtained in a trivial manner 



as 



F LJ r ) = dr'gimw(r,r')sfr c S(r'-r c ) 



C (h^ {ur)ji{ur G )0{T ~ r c ) + 3i{ur)hf- ] (ur c )6(r c - r) 



(D6) 



The radial mode functions induced by the second, third, fifth and sixth terms on the right hand side of (|D5[) are 
obtained in a closed form as 



h m {z) + lA. zh w [z) 



8{r' - r c ) 

\z - \{z - z c )hf\z) + CT^\z)) 0(r - r c ) 

+/ ( ^i(Wc-r)], 

ftT = ~ Z \ (2 3l {zc)h { p{z c ) - i,-i( 2c )C(*) - J l+ i(zc)hf\{z c )) - ^Mzc^izc) , 
" (h{^(zc) - h^zcri^zc)) - \z c hV\z c ) , 



Ml) 

J Imuj 



(D7) 

(D8) 
(D9) 



Fh02(r) = dr'g lmu {r,r')s^l 



0(r' - r c ) 
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1 hm 



(C(z) - hf\ (z)) + \{z- - zl)hf\z) + f^ + h^(z) 9{r - rc) 



M2) + 
J Imuj 

,h(2)- 
J Imu 



Imuj 



(r) = 



Z f (2jK^)M 1) (^)+ij-i(«o)fc&\(^)+im(^)^ ) i(«c)) , 
^{h^\zc) + hf\( Zc )hfUz c )) , 
dr'gi muJ (r, r')sfjnl 

f^l + hi(z)e(r-r c ) 



(D10) 
(Dll) 
(D12) 



i d 

ji(z) + -—zji(z) 
z dz 



6(r c - r') 



•J Imuj 
J Imuj 

F 3{2) (r) 

r lmujV I 



+ [y (3i+i(z) - 31-M) - \{z - zc)ji(z) + ftl~m] 0(rc r) 
Y (j?(zc) - 3i-i{zc)ji+i(zc)) + \zlj}{z c ) , 



(D13) 
(D14) 



(2 J/ (z c )^ (1) (z c ) - ji-^h^zc) - ji + i{zc)hf\{zc)) + ^(z^ (z c ) , (D15) 



= / dr'g lmu (r,r')s^l 



3ii z ) + z—ji(z) 
dz 



0(r c - r') 



„i(2) 



ftTMz)e(r-rc) 



-z (j l+1 (z) + Uz 2 - zl)ji{z) + Ctji(z) } 6{r c - r) 



rJ'(2) + 
J Imuj 

•I Imu 



z: 



f(jf(zc)+3i-i(zc)3i+i(zc)) , 

3 



Z -f {2ji{zc)hf\zc) +ji-x(zc)h\% (z c ) +j l+i (z c )hj 1 ] 1 (zc) 



(D16) 

(D17) 
(D18) 



where we use z — tor and zc — tore- 

Unlike the previous terms, the radial mode functions induced by the fourth and seventh terms on the right hand 
side of (|D5|) cannot be obtained in closed form. For them, it is necessary to obtain the following indefinite integrals, 



£>f )fcfe * = / dzz%(z)kt(z), 



(D19) 



where k and k* represent either j or For finite I, these integrals can be obtained from the reccurence formula 



D 



(l)kk' 



D (ljkk* _ l_ z 2 ( kl{z)kt{z) + kl _ l{z)k! _ l{z) ) 



D 



(2)fefc* 



l + l 
l-l 



-D 



l-l 



1 



2(^-1) 



z 4 (k l (z)kt{z) + k l . 1 (z)kt_ 1 (z)) 



with 



n (l)33 _ 
^1=0 — 


i(ln(z)-Ci(2z)), Djg>i = *- 


z 1 
3 sm(2*)-g 


cos(2z) , 


^1=1 — 


Z —+ Z - sin(2z) + | cos(2z) + ^ ln(z) 


-ici(2z), 




^i=o ~ 


i(Ei(l,-2iz) + ln(z)) , D?2t = 


2 

— + i-e - 

4 4 


1 ilz 

8 ' 


^1=1 — 


£-ife^ + §e^ + ^(1,-24*) 


+ 2^). 




n (l)hh 


Ei(l,- 2i z), D^ h =i\ze^-\ 


e 2lz , 





(D20) 
(D21) 

(D22) 
(D23) 
(D24) 
(D25) 
(D26) 



D 



(2)hh 
1 = 1 



-i-ze 2lz + -e 2lz +Ei(l,-2w). 
2 4 v ' 
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where Ci(z) = — dz cos(z) / z is the cosine integral function and Ei(l, z) — — dze~ zt /t in the first exponential 
integral function. 

Using these integrals, we have 



F^J(r) = dr' 9lmu) (ry)st 3 l 



42 zhf ) {z)+i l ^±^hf\z) 



(r' - r c ) 



1 M3) 



h< L \z)D™ h {z) - Mz)D^ hh (z) - h?\z)Df» h {zc) \9{r - r c ) 
-Mz)Di*> h \zc)6{zc-z) 



F j{3) (r) 



3(3) 
muj 



•n ■ / N -1(1+1) ■ , v" 

-i2z]i{z) + i - J jt{z) 



[rc - r 



J_J(3) 

^2 Imw 



h\ 1 \z)D^(z c )9(z-zc) 



+ hf\z)Df^{z) - ji(z)D^ h (z) + Mz)D^ h (z c ) U(r r c ) 



(D28) 



(D29) 



where we use D\ 



'' : -- -i2D (2)kk * + + l)L> i (1)fc ' c *. 
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